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Abstract 

Computer simulations with Synopsys' Sentaurus TCAD of opto- 
electrical characteristics of separate-confinement heterostructure laser 
based on AlGaAs are used as an example to study the role of the width 
and depth of Quantum Well (QW) active region on laser characteris- 
tics, I-V and I-L, below and above the lasing threshold. The device 
properties depend on both, the number of bound QW states and on 
closeness of the highest bound states to conduction or valence band off- 
set. The lasing action may not exist at certain widths or hights of QW, 
the threshold current is a discontinuous function of these parameters. 
The effects are more pronounced at low temperatures. Discontinuities 
in characteristics may be observed at certain conditions in temperature 
dependencies of laser parameters. 



1 Introduction 



Computer modelling of electronic devices is a relatively new approach to- 
wards study of physical phenomena occurring there as well optimizing their 
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technical characteristics. Methodologicaly, this field of scientific and engi- 
neering activity may be placed between theory and experiment: persuing 
research of that kind requires theoretical understanding of physics of mi- 
croscopic processes and may work as a helpfull tool in interptretation of 
experimental data. However, contrary to sometime met thinking, computer 
modelling can not replace either theory or experiment. But in some situa- 
tions, results of that research may provide an inspiration for understanding 
or testing physical phenomena: it is easier, faster and less expensive to per- 
form modelling than experiments, and we are not restricted that much by, 
often large, inaccuracy of experimental data that may hide insightfull details. 

When performing modelling AlGaAs SCH lasers with Synopsys' Sentau- 
rus TCAD ^Lj, we noticed unexpected steps in some of their characteristics 
(threshold current Ifh versus the width of quantum well da, or versus its 
height, etc). Analyses of results led us to an idea that observed discontinu- 
ities occur when the most upper bound QW state crosses the conduction or 
valance band offset energy. Following that idea, we guessed that the effects 
may manifest itself in temperature dependence of some quantities as well, if 
laser parameters are choosen for that properly. The discontinuities are found 
also below the lasing threshold current, in their I — V characteristics or in 
gain (or loss) versus current. 

We attribute the existence of observed effects to tunneling transport phe- 
nomena through QW. When the dimensions of the structure are reduced 
below the mean free path (which is of the order of 10^ - 10^ nm for AlGaAs), 
one enters the ballistic transport regime; drift-diffusion model of carriers 
transport must be replaced with transfer method type of calculations when 
dealing with QW widths of the order of lOnm [2]. 

2 Modelling 

The laser we are modelling has dimensions, structure and doping as described 
by Andreev, et al. [3j, [4j. The lasing wavelength is 808nm. The lasing offset 
voltage Uq is 1.56 — 1.60V, differential resistance just above the lasing offset, 
r = dU/dl is 50— 80mfi, threshold current Ith is 200— 300mA, slope of optical 
power, S = dL/dl is 1.15 — 1.25W/A, and left and right mirror reflection 
coefficients Ri and Rr are 0.05 and 0.95. The reference laser has the width of 
QW, da, of 12nm and both waveguides' width is 0.2/im. In order to reproduce 
laser characteristics in computational results, we played with several variables 
available in Synopsys. The critical one is Aph - the effective surface area 
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Table 1: A few sets of simulation conditions {A — F) for data shown in 
Figures [2] and m C is radiative recombination rate (Eq. [2]). a„, and ap 
are coefficients of free carrier absorption formula ([T]). Temperature for all 
cases is T = 3007^, electron and hole scattering times are assumed 8 • lO^i^s 
and 4 ■ IQ-^'^s, respectively, and electron and hole mobility 9200cm^/ys and 
AOOcm'^ /Vs, respectively. No SRH recombination and no additional light 
scattering mechanisms are assumed. 



# 


C [cm^/s] 


a„ [cm ^] 


ap [cm ^] 


A 


2.0 • 10-1° 


1 • 10-18 


2 • 10-18 


B 


2.0 • 10-1° 


5 • 10-1^ 


1 • 10-18 


C 


2.0 • 10-1° 


1.5 ■ 10-18 


3 • 10-18 


D 


1.0 • 10-1° 


1.5 ■ 10-18 


3 • 10-18 


E 





1.5 • 10-18 


3 • 10-18 


F 





1.5 ■ 10-18 


3 • 10-18 



factor in Physics section of Synopsys command file. An agreement between 
experiment and calculation is reached for Ap^ of about 0.059. This low value 
of Aph should not be surprising. If carrier drift-diffusion processes only were 
responsible for transport than Aph close to 1 would be expected. However, 
we are dealing here, in particular in waveguide and in QW regions, with 
ballistic transport as well, while the computational model used in Synopsys 
is derived from drift-diffusion equations, modified for dealing with transport 
through QW as discussed in the next section. 

Other parameters available in Synopsys, important in this case, are these 
related to light absorption and carrier scattering. Experiment shows that 
absorption coefficient is of the order of lcm~^ (PI)- It is argued that in 
AlGaAs lasers the main contribution to absorption is due to scattering on 
free carriers. The free carrier absorption coefficient, ajc, is given by: 

a/c = {an ■ n + ap ■ p) ■ L, (1) 
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where n and p are the electron and hole density, and L is light intensity. 
We choosed in our calculations such values of a„ and ap that an effective 
absorption coefficient would be obtained close to that experimental one. 

It is important also to have a reasonable value of radiative recombination 
rate, Rr, which is assumed to be described by: 

Rr = C-(n.p-nl^^y (2) 

where n^^^^ describe the effective intrinsic density, and C is a parameter 
available for changes. 

Typical I — V characteristics computed at T = 300 are shown in Figure [H 
for a broad range of QW widths. For current near the lasing threshold current 
Ith (i-e. for voltage near the lasing offset voltage Uq), which correspond to a 
kink in I — V, for most of these curves the results are very well approximated 
by a phenomenological modified exponential relation ([7J): 

I{U) = Ith ■ exp{A .{U-Uo) + B-{U- Uof), for U <U, 
I{U) = hh ■ exp{C -{U -Uo) + D-{U - U^f), for U > Uo 

where Ith, Uq, as well A, B, C, and D are certain fiting parameters. 

There is no simple, straightforward interpretation of these curves above 
the lasing threshold, where a strong interplay between the effects of carrier 
transport and scattering takes place, with light absorption as well. However, 
we may notice an interesting feature for parts of curves below the lasing 
threshold. While the width of QW, da, changes (nearly) monotonically, the 
curves however are grouped into a few sets such that they nearly coincide 
together within each group. 

A very similar feature is observed when gain or loss is drawn as a function 
of current, for many widths of active region (Figure [3]). Again, we see that 
datacurves for current below the lasing threshold (that corresponds to kinks 
in curves) are grouped into a few sets such that they nearly coincide together 
within each group. If gain or loss were drawn as a function of voltage, 
however, we would not see such a grouping. 

Therefore, we conclude that below the lasing threshold, current as a func- 
tion of QW width at constant voltage derived from the data like these in Fig. 
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[H or gain or loss as a function of QW width, at constant current values, also 
below the lasing threshold, will follow step-like functions. 

This is illustrated in Fig. [2], where current as a function of QW width 
derived at constant voltage from datacurves similar to these as in Figured] is 
shown. Figure [2] presents data computed at different conditions, and marked 
from A to F, for several combinations of free carrier scattering coefficients, 
On and Op, and values of C, the radiative recombination parameter, as de- 
scribed in Table [TJ Moreover, dataset F differs from datasets A-E. The 
last are computed assuming changing Al concentration in waveguides (when 
Al in QW is kept constant) in such a way that the lasing wavelength does 
not change with the change of QW width. The dataset F is computed for 
constant Al concentration in waveguides of 33%. The solid line in Fig. [2] is 
drawn through datapoints F, and arrows there refer to curve F as well, and 
indicate positions of bound QW energy states crossing the conduction- or 
valance band offset energies. Positions of bound QW energy states for cases 
A-E, as illustrated in Fig. [5l are very close to but not identical. 

The step-like features are preserved also in Ith{da) dependencies, as il- 
lustrated in Fig. m There, however, effects of carrier scattering and light 
absorption smear-out the picture. It is useful to notice that at some values of 
QW width no lasing action is reached, and therefore the datapoints in that 
Figure are not available for all QW widths. 

Changes of QW height (caused by difference of Al concentration in QW 
and waveguides) cause very similar step-like dependencies. Moreover, the 
effects are in some situations more clear and pronounced at low temperatures. 
We did modelling for T=77.6 K to confirm their existence. 

Moreover, with a careful design of laser structure (content of Al in QW 
and waveguides) it is possible to find the evidence of the effect in temperature 
dependence of current, when measurements are performed at constant volt- 
age. Energy gap in active region, as well energy gap in waiveguides change 
in a nearly the same way when temperature changes, but not in exactly the 
same way. We found such Al concentrations when the number of QW bound 
states changes with temperature swap. Figure [H] shows how the uppermost 
bound electron state energy, E^^ in this case, differs from Ecbo (energy levels 
of other QW bound states do not play a significant role in this case), for 
a three Al concentrations in waveguide, when Al concentration in QW is 
8%. For Al concentration 34.70%, the E4, energy level exists always through 
temperature swap studied. For Al concentration 34.65%, the E^ energy level 
does not exist between around 150 and 370i^. For Al concentration 34.60%, 
it exists at temperatures higher than about 470i^, only. This has profound 
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implications on /(T) dependencies measured at constant voltage for these 
three different Al concentrations of Al in waveguide, as Figure [7] illustrates. 
For 34.60% and 34.70% of Al content, we observe continuous /(T) curves. 
However, for 34.65% of Al content, at low temperatures /(T) results fall on 
curve that has been computed for 34.60% of Al, and at high temperatures 
they fall on the curve computed for 34.70% of Al. 

3 Discussion 

In case of tunneling energy barrier, transfer matrix approach is used to de- 
scribe charge transport through it ([5j, [6j). The interband tunneling current 
is written as 

J - Ie::: •^(^) • /(^) ■ nE) ■ dE (4) 

where T{E) is energy-dependent tunneling rate, N{E) is the density of 
states, f{E) is the Fermi-Dirac distribution function, respectively, and -Emm 
and Emax are minimum and maximum carrier energies available. 

In Sentaurus, the following intuitive model is used to handle the physics of 
carrier scattering at the quantum well (Figure[8]). The carrier populations are 
separated into bound and continuum states, and separate continuity equa- 
tions are applied to both populations. The QW scattering model accounts 
for the net capture rate, that is, not all of the carriers will be scattered into 
the bound states of the quantum well. The electron capture rate from the 
continuum (subscript C) to the bound (subscript B) states is: 



R = J^ dEc dEB ■ Nc{Ec) ■ Nb{Eb) ■ S{Eb, Ec) ■ fc{Ec){l - fsiEB)) 

(5) 

where E^. and E^ is energy of lowest conduction band-, and bound QW 
electron states, N{E) is the density-of-states, S{Eb,Ec) is the scattering 
probability, and f{E) is the Fermi-Dirac distribution. The reverse process 
gives the electron emission rate from the bound to continuum states: 
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M = J^ dEc dEs ■ Nc{Ec) ■ Nb{Eb) ■ S{Eb, Ec) ■ fB{EB){l - /c(^c)) 

(6) 

The net capture rate is C = R — M, and for very deep quantum wells 
(keyword QWDeep must be used for that in Sentaurus) is known to be given 
by approximation: 

C = R-M = {1- expivB - Vc)) ■ — (7) 

r 

where rjB = {—q^B — Ec) /kBT and rjc = {—q^c — Ec)/kBT contain the 
quasi-Fermi level information and r is the capture time. The capture time 
represents scattering processes attributed to carrier-carrier and carrier-LO 
phonon interactions involving bound quantum well states. The net capture 
rate C is added to the continuity equations as a recombination term. 

In a similar way scattering of holes is computed, with their own char- 
acteristic capture time. These parameters are specified in Sentaurus by the 
keywords QWeScatTime and QWhScatTime. Their default values, 8-10~^^s 
and 4- 10~^^s, respectively, correspond reasonably well to these based on the- 
ory ([8], [9]). Photoluminescence spectroscopy results give values of an order 
of 3 — 20ps [9]. In most of our modelling, if not indicated otherwise, we use 
also default values of electron and hole mobility, represented in Sentaurus by 
parameters eQW Mobility = 9200cm^/Vs and hQW Mobility = 400cm^/V^s. 

For shallow quantum wells, the energy transfer during scattering can only 
occur in a limited range. In the limit of elastic scattering, the net capture 
rate is then approximated by: 

^ ^ ( F3/2ivc) _ F3/2{r]B) \ nc 

where is the Complete Fermi-Dirac integral of the order of m. The 
shallow quantum well model is activated by the keyword QW Shallow. 

It should be pointed out that Equations [7] and [HI for deep and shallow 
quantum wells, respectively, while provide a convenient, intuitive description 
of carriers scattering and capture on QW, these are approximate only. In 
particular, there is no dependence of capture time on energy of unbound 
carriers there and no periodic oscillations as a function of QW size ([8j, [9j). 
GaAs/AlGaAs are considered to have deep quantum wells. However, as we 
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have shown, the resuhs of our modelhng indicate on a strong role of bound 
QW states located very closely to the offset energy levels of quantum wells. 
For these reasons, we did not restrict our calculations to deep- or shallow- 
QW models but used instead the full model available in Sentaurus. 

Equations m - 16] all depend on density of bound states in QW. We expect 
hence that current through the QW will be proportional to the density of all 
bound states in QW. In effective mass approximation, the two-dimensional 
density of electron states within each QW subband n equals ([6]) 



iV„(E) = ^, forE>E^ 

Hence, the current should be proportional to the number of bound states x 
carrier mass. The quantity computed this way (with a certain multiplication 
factor) is represented by large circles in Figure [2l Though it must not be 
exact (for instance, no difference in scattering rates for electrons and holes 
is accounted for), it fits reasonably well I (da) dependence. 



4 Summary and Conclusion 

When performing modeling of laser characteristics as a function of the width 
of active region we noticed a non-monotonic, discontinuous dependence of 
I (da) (when measured at constant voltage applied). A careful analysis of the 
data led us to the hypothesis that discontinuities occur when the most upper 
QW, bound energy states are found very close to the conduction or valence 
band energy offsets. The effect, hence, is thought to be related to changes 
in density of states of carriers from one hand, and to fast changes in carrier 
transfer matrix through QW for QW bound states close to Ecbo oi' Evbo- 
As such, it ought to be more pronounced at lower temperatures. Indeed, 
results of modeling I{da) at liquid Nitrogen temperature (77.6 K) confirmed 
this idea. 

The effect is observed also when modelling current as a function of QW 
depth (Al concentration in waveguide). 

Therefore, we concluded that a similar effect will be present also in mod- 
elling / as a function of temperature. In that case however a carefuU design 
of laser properties is needed, in such a way that a transition of the most 
upper QW energy state will pass through an endge of quantum well when 
temperature is swapt. 
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These observations are important for proper designing of semiconducting 
lasers (choice of Al concentrations, thickness of the active region, etc). How- 
ever, they also illustrate well the intrinsic carrier transport mechanisms in 
SCH lasers. Potentially, might be useful as a kind of quantum level spec- 
troscopy tool when testing laser designs for technological applications. 
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Figure 1: Typical I — V characteristics computed at T = 300. OpticalLoss 
parameter is assumed 0, no radiative recombination, free carrier scattering 
rate parameters r„ and Tp are 8 • 10^^'^s^^and 4 ■ lO^^'^s"^, with electron and 
hole mobolities 9200cm'^/Vs and 400cm^/ys. The legend describes width 
of QW (in nm), from 5nm from right-bottom curve to 24nm for uppermost 
curve. 
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Figure 2: Current as a function of QW width derived at constant voltage 
from datacurves as these shown in Figure [1] and described in Table [T] The 
solid line for datapoints F is to guide the eyes, only. F is computed for 
constant Al concentration in QW of 8%, while all other datasets {A — E) 
are computed with such a concentration of Al in QW that lasing wavelength 
will remain constant (808nm) when QW width changes. The arrows are for 
datapoints F. 
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Figure 4: Lasing threshold current as a function of QW width for datasets 
as these shown in Figure [T] and described in Table [T] The arrows are at 
positions close to but not identical to these in Figure [2J 
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Figure 5: T=300K. Conduction and valence band offset energy, {Ecbo and 
EvBo)i and electron {En)-, light- and heavy hole energies {LHn and HHn) in 
QW, as a function of QW width. Hole energies and Evbo have been scaled 
up by a factor 28 to obtain coincidnce with electron energy scale {Ecbo and 
Evbo curve are the same in this Figure). 
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Figure 6: The energy difference between conduction band offset energy Ecbo 
and the closest to it bound electron state in QW, i?4, as a function of tem- 
perature, when active region concentration is 8% of Al, for three values of 
waveguide Al concentraion: 34.60%, 34.65%, and 34.70%. 
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Figure 7: Temperature dependence of current for free values of voltage ap- 
plied (1.510V, 1.613V, and 1.687V; for each of these this is below the lasing 
threshold), and for free values of Al concentrations in waveguide, the same 
as in Figure |6l 
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Figure 8: The model of carrier scattering at the quantum well used in Sen- 
taurus (based on [1]). 
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